Dephasing due to Intermode Coupling in Superconducting Stripline Resonators 
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The nonlinearity exhibited by the kinetic inductance of a superconducting stripline couples 
stripline resonator modes together in a manner suitable for quantum non-demolition measurement 
of the number of photons in a given resonator mode. Quantum non-demolition measurement is 
accomplished by coherently driving another resonator mode, referred to as the detector mode, and 
measuring its response. We show that the sensitivity of such a detection scheme is directly related 
to the dephasing rate induced by such an intermode coupling. We show that high sensitivity is 
expected when the detector mode is driven into the nonlinear regime and operated close to a point 
where critical slowing down occurs. 

PACS numbers: 42.50.Gy, 42.65.Yj, 42.50.Dv 



I. INTRODUCTION 

The resistive and inductive nonlinearity exhibited 
by superconducting striplines and microstrips originates 
from the nonuniform distribution of the microwave cur- 
rent over the transmission line cross section. Along the 
edges, where the current density obtains its peak value, 
the current density can become overcritical, even with 
relatively moderate power levels. As a result, the super- 
conducting current distribution may change as the total 
current within the transmission line changes. This, in 
turn, causes both the inductance L and resistance R per 
unit length to become current dependent 0. 

Such nonlinearity introduces coupling between dif- 
ferent modes in a stripline resonator. In the ro- 
tating wave approximation (RWA) such coupling adds 
a term to the Hamiltonian of the system given by 
E„y„« ftK'n"N n >N n », where N n = A\A n is the num- 
ber operator of mode n, and the coupling constants \ n 'n" 
are given by IA22I As was shown in Ref . 2] , such inter- 
mode coupling may allow a non-demolition measurement 
> °f the number of photons in a signal mode by in- 
tensively driving another mode, called a detector mode 
(or pump), and monitoring its response near the pump 
frequency. Such a measurement scheme is character- 
ized by a measurement time, defined as the time needed 
to distinguish between initial states of the signal mode 
having different number of photons. Single-photon de- 
tection can be realized if the measurement time can be 
made shorter than the lifetime of a Fock state of the sig- 
nal mode. 

However, according to the Bohr's complementarity 
principle [5|, in the limit of single-photon detection, de- 
phasing is expected to come into play, leading to broad- 
ening of the resonance line shape of the signal mode. 
Thus, a which- path-\\ke experiment can be employed to 
test whether single-photon detection is possible. This 
can be done by monitoring the resonance response of the 
signal mode and the way it is affected by driving the de- 



tector mode. Observing a broadening comparable to or 
larger than the width of the resonance when the detec- 
tor is not driven implies that single-photon detection is 
possible in principle, namely, the measurement time is 
comparable to or shorter than the lifetime of a photon in 
the mode. 

In practice, however, as we show below, the coupling 
constants \ n i n " in superconducting stripline resonators 
are far too weak to allow single-photon detection when 
the detector mode is taken to have a linear response. 
On the other hand, the same kinetic inductance which 
leads to intermode coupling also gives rise to a detector 
mode Kerr nonlinearity. In the RWA such nonlinear- 
ity adds the term J2 n K n (N n ) 2 to the Hamiltonian of 
the system, where the Kerr constants K n are given by 
IA23I Such Duffing-likc nonlinearity leads to bistability 
when the drive level exceeds some critical value. In the 
bistable regime of operation the response of the detector 
mode exhibits hysteresis and jumps as the frequency of 
the drive or its amplitude is varied. 

In the present work we investigate the dephasing rate 
induced by the driven detector mode on the signal mode. 
We consider a detector mode having both a Kerr-likc 
nonlinearity and nonlinear damping. We find that the 
dephasing rate diverges as one approaches a point where 
the response of the detector mode exhibits a jump. At 
such a point the slope of the response vs. frequency is in- 
finite. This result suggests that strong dephasing can be 
induced when driving the detector mode in the nonlinear 
regime, even with relatively weak intermode coupling. 
Thus, operating such a photon detector in this regime 
may allow enhanced sensitivity. 

The dynamics of a driven single-mode stripline res- 
onator having a Kerr-like nonlinearity, in addition to 
linear and nonlinear damping was studied recently 
by us in 0. The equations of motion were derived 
using the input-output theory of Gardiner and Collett 
0, U We have shown that in the appropriate limit such 
a system can serve as a phase-sensitive amplifier whose 
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noise performance exceeds that of the quantum limits im- 
posed on linear, phase-insensitive parametric amplifiers 
[lCj . Moreover, we have studied degradation of device 
performance, due to two-photon absorption originated 
by nonlinear resistance, which can be significant in such 
devices (see also [H El ES El ) . 

A closely related problem, the non-demolition mea- 
surement of Fock states in mesoscopic mechanical oscil- 
lators, was recently studied theoretically in Anhar- 
monic effects on such a measurement were further studied 
in [l^. In the RWA the Hamiltonian of the mechani- 
cal systems, studied in E3 an d 0i an d the one 01 a 
stripline resonator presented here, are similar. However, 
our analysis goes beyond that of 15] and fill bv taking 
into account the effect of nonlinear dampin g |l7j . More- 
over, while the emphasis in Refs. |15| and [16j is on the 
properties of the detection scheme, including a thorough 
consideration of the phase diffusion back-action, here we 
focus on the relation between decoherence and the dis- 
tinguishability between states with different numbers of 
photons. 



II. SINGLE-PHOTON DETECTOR 

The Hamiltonian of the system is given by 

H = H S +H D + V . (1) 

Here H s is the Hamiltonian of the mode of the photons 
to be detected (signal mode), 



(2) 



where N s — A\A S is the number operator of the signal 
mode. Hd is the Hamiltonian of the detector. V is the 
coupling between the signal mode and the detector, and 
is taken to have the form 



V = h\N s W, 



(3) 



where W is an operator on the Hilbert space of the de- 
tector. 

The Heisenberg equation of motion can be used to 
show that the number of photons in the signal mode N s 
is a constant of motion 



\N S ,H] 



0. 



(4) 



Thus, the detection process in this case is a quantum 
non-demolition measurement of the number N s • 

To study the ability of the detector to measure N. w e 
follow the approach presented in Refs. E3i |20f. 
Consider an initial state at time t = 0, having a super- 
position of and 1 photons 



l^(* = 0)) = (|0), + |l),)(8|xi> 



D I 



(5) 



where \n) s is a Fock state of the signal mode and \xi) d 
represents an initial state of the detector. At a later 



time t > the state in the Schrodinger representation 
will evolve into 

|^(t))=tt(t)|^(t = 0)), (6) 
where the time evolution operator u (t) is given by 



(t)=T exp 



dt'H (t r ) 



(7) 



and T is the time-ordering operator. Equation ^ yields 

u(t)\n) s ®\xi) D (8) 



= Texp 



-~ / dt [H D + h(u> s + XW) N s 
ft Jo 



\ n )s ® \Xi)l 



= \n) s ®u n \Xi) D , 
where 



— inojatrr" 

u n = e 3 i exp 



D + nhXW) 



(9) 



Thus, as expected, the time evolution preserves _/V s , and 



hM*)> = |o>, ® Ixoh 



i 



IXi)d, 



(10) 



where \xo) D = u \xi) D and \xi) D = ui \Xi) D - The dis- 
tinguishability between and 1 photons associated with 
this state is characterized by the parameter 



v = \(A s )\' 



(Xolx 



D 



(X*l4 u i \Xi) D 



(11) 

Our goal is to find the characteristic time scale t v over 
which v decays from its initial value v — 1 at time t = 0. 
Regarding the coupling parameter A as small, one can 
evaluate v using perturbation theory. To lowest order in 
A one finds using IbToI 



v=l-\ 2 Jdt'Jdt"K(t' 1 t"). 
o o 



(12) 



where 



k {t 1 , t") = i [(w (*') w (*")) + (w (t") w (t') 



and 



W(t) =W(t)- (W(t)) 



(13) 



(14) 



For a steady state the correlation function K(t',t") is a 
function of r = \t' — Moreover, it decays on some 
characteristic time scale t v , which can be identified from 
Eq. [H 19] 



— = A" 



/oo 
drK (t) 
-OO 



(15) 



where K (t) =K (0,t). 
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III. THE DETECTOR MODE 

The detector mode is a resonator mode distinct from 
the signal mode. Because it is driven to large amplitude, 
the Kerr nonlinearity and the linear and nonlinear damp- 
ing of the mode are taken into account in our analysis. 
Here we review the main results of Ref . . 



A. The Hamiltonian 

The Hamiltonian is given by 

Hd = Ha + Hk + H ai + H ao + H a3 + Ht x + Ht + Ht 3 , 

(16) 

where Ha is the Hamiltonian for the detector mode 

H A = hw^A , (17) 
Hk is the Hamiltonian for the Kerr nonlinearity 



H K = -KA^A^AA , 
2 



(18) 



H a \, H a 2 and 7J q3 are the Hamiltonians for the bath 
modes associated with the dissipative elements 



H a i = / dojhwa\(uj)ai(Lu) 



(19) 



B. Equations of motion for A 

The equation of motion for A is given by 

— = -iuj Q A - iKA^AA --/A- j 3 A^AA (25) 
at 

- i2^e^AUf(t), 

where 7 = 71 + 72 . The damping constants are written 
as 



where j n and <p n [n = 1, 2, 3) are real, and 
1 



(26) 



<{t) 



doje- w{t ~ ta) a n {t ,uj). (27) 



The following commutation relations exist for the bath 
modes 



[a n (t ,uj),al n {t ,uj')\ =S(u - w')S v 
[a n (t ,uj),a m (t ,uj')} = , 



(28) 
(29) 



where n, m = 1, 2, 3. 



Mean-Field Solution 



H a 2 = J <kufojja\{uj)a,2{<jj) 



H a 3 = / du!hxjja\{ui)a 3 (u>) 



(20) 



(21) 



Ht x and Ht 2 are the Hamiltonians for the linear dissi- 
pative elements, one of which is the port through which 
signals enter and leave the oscillator. These Hamiltonians 
linearly couple the bath modes, a\ and 02 respectively, 
to the oscillator mode A 

H Tl =h J duj[TxA^ a x {u) +T*a\{u)A} , (22) 



H T2 =hJ <kj[T 2 Aio2(w)+T£al(u))A] 



(23) 



The two-photon absorptive coupling of the resonator 
mode to the bath modes 03 is modeled by a hopping 
Hamiltonian in which two cav ity pho tons are destroyed 
for every bath photon created [21II2II2II2II2H 



H T3 =hJ du{T 3 A^A^a 3 (cj)+T^al(cj)AA}. 



(24) 



The modes are boson modes, satisfying the usual Bosc 
commutation relations. 



Since the detector mode is driven to large amplitude, 
compared with quantum fluctuations, the behavior of the 
mean-field is the classical behavior of the mode, except 
very near instability points. Hence, to obtain expressions 
for the mean-field behavior, in this section, we treat of 1 , 
a™, and a™ as complex numbers rather than as opera- 
tors. We take the bath mode amplitudes a\ n and a 3 ™ to 
be those entering the cavity from losses. Consequently, 
we take a\ n = and a 3 ™ = 0. The amplitude af 1 is cho- 
sen to be that of the incoming drive and is taken to have 
the oscillatory time dependence 



,in _ ijn -i(u> p t+4>i) 



(30) 



where b\ n is a real constant. Writing A as 

A = Be^ Wpt+ ^ B \ (31) 

where B is a positive real constant, the equations of mo- 
tion become 

[i(w -w p )+7] J B+(iA>7 3 )B 3 = -iV^re^ 1 ^-^ 1 ' ■ 

Multiplying each side of the nonlinear equation by its 
complex conjugate, one obtains |2(j, p7| 



2[(w - uj p )K + 773] 



K 2 - 

{ujq - uj p ) 2 + 7 2 2 

H „o T & — 



7i 



271 



B 4 



(33) 



K 2 



7l 
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Once B 2 has been determined from the above cubic al- 
gebraic equation, the phase can be found from 

oot(fr + 4> B - V>i ~ tt/2) = 7 + ^ R2 ■ (34) 

Taking the derivative of Eg. 1331 with respect to u p , one 
finds 



dB 2 
duj n 



(35) 



2(cj - uj p + KB 2 )B 2 



(lo q -lo p + 2KB 2 ) 2 + ( 7 + 2 l3 B 2 Y - (X 2 + ~/ 2 )B* 



where £? constitutes the mean-field amplitude of the de- 
tector mode in response to the classical drive b™ . The 
operators c"\ c™, c|™, and a are regarded as small and 
will be kept only up to linear order. In this approxima- 
tion the operator a satisfies the following second-order 
equation of motion 

^ + 2R(«,)^ + (H a -H 2 )a = r(t), (45) 



where 



w = i(Lj -uJp)+j + 2{iK + j 3 )B 2 , (46) 



D. Onset of testability point 



V = (iK + j 3 )B 2 e 



2„-2iA 



(47) 



At the onset of bistability point the following holds 



dui p _ d 2 uj p _ 
dB 2 ~ d(B 2 ) 2 ~ 



(36) 



Such a point occurs only if the nonlinear damping is suf- 
ficiently small 



\K\ > x/3 73 . 
At this critical point the following holds 

B 2 - 27 
c " V3(\K\ - a/373) ' 



(37) 



(38) 



T(t) = ^-+w*F(t)-vF^t), (48) 

and 

F (t) = -iv^e^c™ - iV272e^ 2 4" (49) 
The solution to Eq. 03] is given by 

/oo 
drG(t - r)r(r) , (50) 
-00 

where 



^0 - ^pc = -7- 



A 
^1 



473|A| + V3(A 2 + 73 2 ) 
A 2 -37? 



{b in ) 2_ 4 7 3 (A 2 +7l) 



3x/3 7 i(|A|- V3 73 ) 3 



(39) 



(40) 



G(t) = u{t) 



Q — ^ot g— Ait 



Aq — Ai 



(51) 



u(t) is the step function, and the eigenvalues Ao and Ai 
satisfy 



A + Ai = 2S(u>), 



(52) 



E. Quantum fluctuations about the mean-field 
solution 



To determine the quantum fluctuations about the 
mean-field solution we write 



and 



in _ fin —i(Ljpt+il>i) 1 in 



— iuj p t 



in _ in -iuipt 
"2 — c 2 c i 



in „in„— iw D t 



a 3 = c 3 e 



A = Be- l{u> " t+4 ' B) + ae^K 



(41) 
(42) 

(43) 
(44) 



A Ai = \w 



2 -hi 2 



(53) 



From Eq. [23J one finds that the slope of the response 
function B 2 vs. ui v is infinite when 



\w\ 2 -\v\ 2 = 0. 



(54) 



Thus, using Eq. [331 one finds that at these points at least 
one of the eigenvalues Ao and Ai vanishes. The system 
eigenmode corresponding to the eigenvalue that vanishes 
is the one that experiences critical slowing down. 



F. Thermal noise 

Assuming the bath modes are in thermal equilibrium, 
the following hold 



<C(r)> = (C t (r)>=0, 



(55) 
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(C(t)C(t')) = (CHt)CHt')) = 0, (56) 

(4 nt (r)C(r')) = * (r - r') £ nm (n WD ) , (57) 

(C(t)CV)) = J (t - r') «5 nro «n Wo ) + 1) , (58) 
where 

= ^ri' ( 59 ) 

n, m = 1,2,3 and /3 = l/ksT. 

IV. THE DEPHASING RATE 

Since the interaction Hamiltonian between modes is 
bilinear in the number operators for the modes, one has 
W = A^A, or using Ea.ltll 

W (i) = B 2 + B [a (t) e^ B + a f (t) e"^ B ] + a f (i) a (t) . 

(60) 

Using Eas.l5UI and 1551 one finds 



(a(i)) = (« f (t)) = 0. 



(61) 



Moreover, keeping terms only up to first order in a, one 
finds 

W (t) = W (t) - (W (t)) =B[a (t) e l(t>B + a f (*) e^ B ] . 

(62) 

In order to use Eq. El to calculate the dephasing rate, 
some expectation values must be evaluated. Using 
Eqs.EHlEniEZI andEHl one finds 



(F(r)F(r'))=0, 



(Ft (r)Ft (r')) = 0, 



(63) 



(64) 



(F(r)Ft ( T ')) = (Xo+Xi)S(T-r')(n Uo ), 



and 



(Ft (r) F (r')) = (A + Ai) 5 (r - r') ((^ ) + 1) . (65) 
Thus, using Eq. 023 one finds 
(a(t)a(t')) (66) 

pOO 

dr I dr'G(t - r)G(t' - r') (T(t)T(t')) 
dG(t-r) 



— OO <J oo 



= (A + Ai) (n Uo ) v dr 



dr 



-G(t'-r) 



+ (A + Ai) ((n Wo ) + l)v / drG(t - r) 



( dG{t' - t) 
~d~r 



(A + Ai) (2 (n U0 ) + 1) uw* / drG(t - r)G{t' - r) 



and 

( O t(t)o(0> (67) 
//r / dr , G(t-T)G(t'-r , )(rt(T)r(r')) 

dG(t - r) <fG(t' - t) 



= (A + Ai)((n wo ) + 1) 



dr- 



dr 



dr dr 
dG(t - t) 



dr 



/OO 
drG(t - t) 
-OO 



G(t'-r) 
dG(t' - t) 

+ (A + Ai) ((n^) + 1) \w\ 2 I drG(t - r)G(t' - r) 

J — OO 

/OO 
drG(t - r)G(t' - t). 
-OO 

Using Eq. these results yield 
K (i, t') 



B 2 (A + Ai)coth 



(68) 



dr 



dG{t - t) dG(t' - t) 



Re (w — ve 



dr 

2i<j> B 



dr 



G(t 



dr 



dGlt' - t) dG(t - t) 



dr 



dr 



G(t' 



+ \ w — ve 



2i<ji E 



' / drG{t - T )G(t' - t) 

J — OO 



From this, by noting that, from Eq s . 1461 ET7I and l3*2l 



(w - ve 2 ^ B ) B = -U 
and that 



i(<t>i+4>B-i>i) 



(69) 



/OO pOO -i 

dt' / drG(-r)G(t' - r) = —— 2 
-oo J-oo (AoAl) 

one obtains, from Eq. 1151 the dephasing rate 
-j 2 7l (6 X ) coth— — . 



1 _ ^2 Aq + A i 
T<P (AqAi; 



(70) 



Thus, 1/t v is directly related to the eigenvalues Ao and 
Ai which characterizing the response of the detector 
mode to small perturbation. As was shown above, 
the product XqXi vanishes at the points where the slope 
dB 2 / dujp is infinite, leading to a diverging dephasing rate 
1/t v . However, our model, which takes nonlinearity into 
account only to lowest order, breaks down near these 
points, where the fluctuation in ^4t^4 becomes apprecia- 
ble. To evaluate the actual dephasing rate near these 
points one has to take into account higher-order terms 
in the nonlinear expansion. This is beyond the scope of 
the present work, although we note that the divergence 
of X/tu, in the present model indicates that the actual 
dephasing rate will be relatively large near these points. 
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FIG. 1: Square root of average photon number in the de- 
tector mode B and log of the normalized dephasing rate 
log (7/A 2 t^) , for subcritical drive b 1 " = 0.5&i™ (panels (a) 
and (b)), critical drive 6| n = &!" (panels (c) and (d)), and 
overcritical drive = 26™ (panels (e) and (f)). The pa- 



rameters in this example are K 
72 = I.I71, and 73 = W~ 2 K/V3. 



-10 u>o, 71 = 10 wo, 



For the linear case K — 0, 73 = the dephasing rate 
is given by 



1_ = 4A 2 77l (61") 2 
Tip [(cj - uj p ) 2 + 7 2 



■ coth 



(71) 



Figure ^ shows B, found from solving Eq. and the 
log of the normalized dephasing rate, log (j/X 2 ^), cal- 
culated by using Eq. [7U1 Three cases are shown, sub- 
critical drive 6™ = 0.56™ (panels (a) and (b)), critical 
drive 6 m = 6™ (panels (c) and (d)), and overcritical drive 
b" 1 = 26™ (panels (e) and (f)). The values chosen for 
this example for the parameters K, 71, 72, and 73 are 
indicated in the figure caption. As can be seen from 
Fig. n (d) , the dephasing rate is an asymmetric function 
of frequency around the onset of bistability point, where 
1/t v diverges. Also, comparing the two jump points in 
Fig. ^ (f) indicates that the divergence in the right one 
(with higher frequency) is steeper. 



V. DISCUSSION 

The sensitivity of single-photon detection can be char- 
acterized by the dimensionless parameter £, defined as 



1 

7V 



(72) 



Consider first the linear case K = 0, 73 = 0, where the 
dephasing rate is given by Eq. 1711 We assume for sim- 



plicity 271 ~ 7, and consider the case of low temperature 
(3Hujq 3> 1. At resonance where the dephasing rate is 
largest, one finds 



< = -IB- ( ^ 



(73) 



In the present case of linear response we assume B < 
B Cl where B c is given by Eq. |2H1 Moreover, as can be 
seen from Eq. (|A22() and l|A23(l . in many cases K is of 
the same order as A. Thus, assuming also the case where 
73 « K, one finds that the largest possible value of £ 
in this regime is roughly given by 



A 
7' 



(74) 



The experimental data in Refs. [2|| and |22|, together 
with Eq. |2H1 allow one to estimate the value of K. As- 
suming, as before, that A ~ \K\, one finds for striplinc 
resonators made of Nb that £ max — 10~ 12 . For res- 
onators made of NbN, on the other hand, the onset of 
nonlinear instability occurs at a much lower input power. 
However, the nature of the nonlinearity observed in NbN 
resonators is qualitatively different from what is expected 
from nonlinearity originated by the kinetic inductance 
effect, thus indicating that the underlying physics is dif- 
ferent |29| . If, however, in spite of this discrepancy, we 
employ Eq. [JElalso for the case of NbN resonators, we ob- 
tain £ max ~ 10~ 6 . We thus find that for both cases £ max 
is far too small to allow single-photon detection when the 
pump mode is operated in the regime of linear response. 
In contrast, as was mentioned above, operating the pump 
mode in the regime of nonlinear response may allow a sig- 
nificant enhancement in the sensitivity. Further study is 
required, however, to analyze the behavior of the system 
close to the points where our model yields divergence in 

i/v 



APPENDIX A: LOSSLESS TRANSMISSION LINE 
RESONATOR 

Consider a lossless linear transmission line with length 
I extending along the x-axis 0. Let q (x, t) be the charge 
density per unit length and define 



Q(x,t) = I dx'q(x',t). 



(Al) 



Thus, q = —dQ/dx and the voltage across the transmis- 
sion line is given by 



(A2) 



where C is the capacitance per unit length along the 
transmission line. The current is given by 



I{x,t) 



dQ 
dt ' 



(A3) 
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The Lagrangian of the system is given by 



and 



C= \ I dX ~ CV ^ 



(A4) 



I 



1 



= 2 ldX 

o 



5Q\ 2 1 ( dQ s 1 
~dt 



C \dx 



where L is the inductance per unit length along the trans- 
mission line. Here we assume that both C and L depend 
on x. Moreover, the kinetic inductance leads to a de- 
pendence of L on the current, given by pj, 



L = L + AL 



I 



As a basis for expanding Q (x, t) as 

( x > t) = ^ <7n (*) «« («) , 



(A5) 



(A6) 



we use the solutions of the boundary value problem 

^(c^)=-^« (A7) 

with the boundary conditions of vanishing current 
u n (0) = u n (I) — 0. We assume that the functions 
u n (x) are chosen to be real. Using this basis and the 
Strum - Liouville theorem, one obtains 



C = \ E (9n " "111) + A£ ' 



where 



n' ,n ,f ,n /,f ,n /,f 



Qn' Q_n" yln'" Qn 



(A8) 



(A9) 



X / dxALUn'U n iiUn'"U n "" ■ 




The variable canonically conjugate to <7„ is 



<9£ 



9A£ 



To first order in AL the Hamiltonian is given by 



(A10) 



where 



Ho = jE(^+^) (AH) 



V = -^J2 Pn'Pn"Pn"'Pn"" (A12) 



n' ,n" ,n'" .n" 



x / dxALu n 'U n "U n '"U n "ii . 



Quantization is achieved by regarding the variables 
{gVuPn} as operators satisfying the following commuta- 
tion relations 



[?niPm] = q n pm - p m q n = ihS n 



(A13) 



and 



[q n ,q m ] = \Pn,Pm] = . (A14) 
In terms of the Boson annihilation and creation operators 



A n = ( V^VM + ~l=Pn 



„ . „ W„9„ —p n , 

2ft V ^^^^ 



(A15) 



(A16) 



the Hamiltonian Eq. I|A11|I can be expressed as 

Ho = ( A l A n + I) ■ (AIT) 



The current operator is given by 
I(x,t) = ^=i 



^2\l^(Aie^-A n e-^) Un (x) 

(A18) 



and the voltage operator is given by 



n 

(A19) 

The factors p n 'Pn"Pn"'Pn"" generally contain terms os- 
cillating rapidly at frequencies on the order of the fre- 
quencies in the resonator spectrum. In the RWA these 
terms are neglected, since their effect on the dynamics 
on a time scale much longer than a typical oscillation pe- 
riod is negligibly small, and only stationary terms remain 
(lfij . Thus, in the expression of V only terms of the type 
Pn'Pn" contain stationary terms which are given by 



2 2 hu) n i fth) n " 
Pn'Pn" 



1 + 2A\ l ,A n ^j (l + 2A\,,A n , 



2 2 

(A20) 

The constant term can be disregarded, since it only gives 
rise to a constant phase factor. Moreover, the terms 
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X* n ,A n i and A\,,A n n that give rise to frequency renor- 
malization can be absorbed into H.q. Thus, in the RWA 
the perturbation V is given by 



V= J2 h\ n , n „Ai,A n ,Ai„A n „ +J2 hK n> (Al,A n ,y , 

(A21) 



n'^n" 

where 



and 



K'n" = --^hj n 'U) n >> J dxALu^u^f (A22) 



K n > = ~^p^l> J dxALu*,. (A23) 



APPENDIX B: PERTURBATION THEORY 



Consider the Schrodinger equation 

i±\i/>) = K\1>), 

where JC = K,^ is given by 

/C = /Cq + A/Ci , 



(Bl) 



(B2) 



where A << 1 is real. The time evolution operator 
u (t, to) can be expanded as |3f)j | 

u (t, t ) = u (t, t ) + Xui (t, t ) + X 2 u 2 {t, t ) + O (A 3 ) , 

(B3) 

where 



Using this expansion, one can calculate the operator 
O (t) = i*q (t, t ) U (t, to) to second order in A 



0(t) 



(B6) 



1 - i\f dt'uo (t , t') JCi (t 1 ) uo (*', to) 



t t' 

X 2 fdt'Jdt"u {to, t') Ki (0 mo (f, t") ACi (t") uo (t", t ) 
to to 



or 

t t t' 

O (t) = l-iXfdt'Ki H {t')-X 2 Jdt'J dt"KiH (f) fC 1H (*") , 

to tQ to 

(B7) 

where 

/C m (t) = 4 (t, t ) ACi (t) M (t, t ) . (B8) 
Since /Ci (t) is Hermitian, one finds to lowest order in A 



\{0(tw 



(B9) 



1 - X 2 Jdtjdt" [{K t (t')/Ci (t")> - (iCi (t')> (ACi (t")>] 

to to 



or 



| (O (t)) | 2 = 1 - X 2 Jdt'Jdt" (Ki (f) Ki («")) , (BIO) 



where 



t x (t) = ACi (t) - (ACi (t)) 



(Bll) 



ui (t, t ) = -»/ dt'uo (t, t') /Ci (t') u (f, t ) (B4) 

to 



and 



u 2 (t,t ) (B5) 
t t' 

= -Jdt'Jdt"u (t, t') ACi (t') u (f, t") /Ci (t") u (t", to) • 

*0 *0 
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